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The phase behaviour of colloidal dispersions is interesting for fundamental reasons and for technological 
applications such as photonic crystals and electronic paper. Sedimentation, which in everyday life is relevant 
from blood analysis to the shelf life of paint, is a means to determine phase boundaries by observing distinct 
layers in samples that are in sedimentation-diffusion equilibrium. However, disentangling the effects due to 
interparticle interactions, which generate the bulk phase diagram, from those due to gravity is a complex 
task. Here we show that a line in the space of chemical potentials ju b where i labels the species, represents a 
sedimented sample and that each crossing of this sedimentation path with a binodal generates an interface 
under gravity. Complex phase stacks can result, such as the sandwich of a floating nematic layer between top 
and bottom isotropic phases that we observed in a mixture of silica spheres and gibbsite platelets. 
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The analysis of the effect of gravity on colloidal dispersions dates back to Perrin 1 and is used in computer 
simulations 2 ' 3 and experiment with e.g. depolarized light scattering 4 " 6 to obtain the osmotic equation of state 
over a wide range of densities from a single sample. Added depletion agents, e.g. non-adsorbing polymers 7 , 
modify the effective interactions between particles of the primary component, but are typically gravity-neutral. 
However, in colloidal mixtures all components are subject to gravity, and hence compete to minimize their 
gravitational energy. The strength of gravity can be quantified by the gravitational lengths ^ = k B T/(mjg), where 
k B is the Boltzmann constant, T is absolute temperature, m f is the buoyant mass of species i (obtained by 
subtracting the solvent background), and g is the acceleration due to gravity. Colloidal liquid crystals 8-11 constitute 
ideal candidates for the investigation of gravitational effects due to their rich phase behaviour. 

We used electro -sterically stabilized mixtures of gibbsite platelets and alumina-coated silica spheres (Klebosol 
30CAL25 and 30CAL50). Particle dimensions were obtained from transmission electron microscopy and atomic 
force microscopy (AFM). The average bare diameter of the spheres was cr s = 30 nm and 74 nm for Klebosol 
30CAL25 and 30CAL50, respectively, with polydispersity of ca. 15% and 21% (3 volume % of the particles are 
below 60 nm in the case of 30CAL50). The average diameter of the platelets was 186 nm with 29% polydispersity. 
AFM measurements of a diluted platelet sample deposited on a mica substrate gave a bare platelet thickness of d ~ 
5 nm. As such results can easily be affected by differences in chemistry between substrate and platelets, we rather 
treat d as an adjustable parameter, with d = 3.7 nm giving the best agreement of the isotropic-nematic (IN) 
coexistence densities obtained in theory 12 and experiment 13 . The value of d enters the conversion from packing 
fraction, r\ b to number density, p b via p { = y\Jv„ where v { is the bare particle volume of species i = S (spheres), P 
(platelets). Both species were purified by dialysis against deionised water containing 5 mM NaCl in order to 
screen electrostatic interactions. Both species had a stabilizer (Solsperse 41000) adsorbed onto them. As a result, 
small-angle neutron scattering showed effective particle diameters o* that were larger by ~ 10 nm than the bare 
sizes, hence a* = a,- + 10 nm. The sample preparation follows that in a closely related system 13 . The (bare) packing 
fraction of spheres was rj s = 0.05 throughout, and we considered (bare) platelet packing fractions rj P = 0.01, 0.025 
and 0.05. The gravitational lengths were £ P = 2.92 mm for the platelets, £ s = 22.4 mm for the small spheres, and 
£ s = 1.49 mm for the large spheres. These were obtained by considering the mass density for silica 2.30 g/cm 3 , for 
gibbsite 2.42 g/cm 3 , and for the aqueous solvent 1.00 g/cm 3 . The density of the adsorbed layer of stabilizer is 
1 .02 g/cm 3 . As the latter value is very close to the solvent density, we neglected its effect on the buoyant masses and 
hence the gravitational lengths of both species. Experiments were carried out at room temperature, T = 293 K. 
The value of the temperature is irrelevant for the slope of the sedimentation path, as this depends only the ratio of 
the gravitational lengths, such that the temperature dependence of the f cancels out. 



SCIENTIFIC REPORTS | 2 : 789 | DOI: 1 0.1 038/srep00789 



1 




Figure 1 | Photographs taken with white-light illumination of samples between crossed polarizers, indicating the time evolution of sedimentation in 
colloidal mixtures of gibbsite platelets with diameter cr P = 186 nm and silica spheres with diameter <7 S = 30 nm (top row) and 74 nm (bottom row). 

The initial bare packing fraction of spheres was rj s = 0.05, and of platelets rj P = 0.01 (left), 0.025 (middle), and 0.05 (right, as indicated). Both sets of three 
samples are shown in the initial state immediately after filling (first column), after 12 hours (second column), after 5 months (third column), and after 20 
months (fourth column) under gravity. The internal width of the vials was 10 mm, the height 40 mm, and the thickness 1 mm. The position of the 
meniscus between the dispersion and the air is highlighted (white lines). The blurred arcs in the lower part of the pictures for rj s = 0.01 and 0.025 (bottom 
row) are artifacts due to the photographic setup. 



Results 

Photographs of the time evolution of the different mixtures are pre- 
sented in Fig. 1 . The initial filling process effectively homogenized the 
systems, with nearly uniform birefringence indicating that nematic 
order formed quickly at short times. After 12 hours the mixtures 
with small spheres had separated into a lower nematic and a super- 
natant isotropic phase. Subsequently, the interface moved down- 
wards compressing the nematic phase. The time evolution of the 
mixtures with the larger spheres was surprisingly different. While 
the sample with y\ p = 0.01 remained nearly isotropic after 12 hours, 
that with r\p = 0.025 had developed a weak gradient in birefringence. 
Only the sample with r\ P = 0.05 showed clear phase separation, albeit 
in reversed order, with a top nematic and bottom isotropic layer. 
After 5 months its IN interface had moved to smaller heights. The 
two other samples, however, had developed a floating nematic layer 
between isotropic bottom and top layers. This state remained stable 
and could be observed after 20 months. Even the sample with r\ P = 
0.05 had transformed into a thick floating nematic. 

In order to rationalize these observations, we use a density func- 
tional theory for a mixture of hard spheres and vanishingly thin hard 
platelets 14 with (effective) diameters cr*. The (Helmholtz) free energy 
functional, F[p s , p P ], reduces to the Rosenfeld functional 15 for pure 
hard spheres, and includes contributions up to third-order in platelet 
density; here p { is the one-body density distribution of species i. The 
theory goes beyond Onsager's treatment 16 , and accurately describes 
the IN phase transition in pure hard platelets 12 . Numerically solving 
the Euler- Lagrange equations, SF[p s , p P ]/dpi = fa, gives the equilib- 
rium densities p- q (^ s ,^ P ), for i = S, P. 

We use the common splitting of the free energy functional into 
ideal and excess contributions, F = F id + F exc . Here the ideal free 
energy functional for the current binary system is F ic i = j drp s (r) 
\\n(p s (r)A 3 s ) — l] + \dx Jdcop P (r,co) [ln(p P (r,co)Ap) — l], where 
r is the position coordinate and w is a unit vector normal to the 
platelet surface that describes the platelet orientation; the spatial 
integral is over the system volume, and the orientation integral is 
over the unit sphere; A/ is the (irrelevant) de Broglie wavelength of 



species i = S,D. The excess free energy functional F exc involves con- 
volutions of both bare density fields with appropriate weight func- 
tions. The resulting weighted densities are combined into a free 
energy density per volume, and integrating this over the system 
volume yields the excess free energy. For further details of F exc we 
refer the reader to the original publication devoted to a ternary 
sphere-needle-platelet mixture 14 and to the recent description of 
the platelet- sphere subsystem 17 , where the binary platelet-sphere 
functional is described in detail. 

Equating pressure, T, and \i { in both phases gives the bulk phase 
diagram, which we show in the density representation in Fig. 2a (b) 
for the mixture with smaller (larger) spheres. The topology of the 
phase diagram is the same in both systems. Pure hard platelets have a 
very weak first-order phase transition that arises from the competi- 
tion of rotational and translational contributions to the free energy: 
in the nematic phase the loss of orientational entropy is over- com- 
pensated by the more efficient packing of the aligned platelets 18,19 . 
Adding spheres widens the density gap considerably, and these pre- 
fer the isotropic phase over the nematic phase. For high densities an 
isotropic state, almost pure in spheres, coexists with a nematic phase 
that is almost pure in platelets. 

We expect that equilibration on small length scales had occurred 
in experiment after 12 hours, and hence represent the average densi- 
ties as statepoints. For the statepoints inside the biphasic region, the 
relative position on the corresponding tieline directly gives the rela- 
tive volumes of the two phases, cf. the cartoons in Fig. 2a and b. The 
top nematic phase in Fig. 2b is a simple consequence of its smaller (in 
this case) total mass density, m s p s + m P p P . Besides the thin nematic 
bottom layer for the small spheres at r\ P = 0.01 and the birefringence 
gradient in the large sphere system for rj P = 0.025, the agreement of 
theoretical and experimental results is remarkable. 

In order to analyze the late stages of equilibration, we turn to the 
chemical potential representation of the phase diagram, cf. Fig. 2c (d) 
for the mixture with small (large) spheres. The biphasic region col- 
lapses to a binodal line. The binodals of both systems are again very 
similar to each other, with a vertical slope for pure platelets, /n s — > 
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Figure 2 | Theoretical bulk phase diagrams of mixtures of hard spheres and infinitely thin hard platelets with aspect ratio 0 P /cr* s = 4.9 (a, c) and 2.33 (b, 
d), corresponding to the experimental mixtures with cr s = 30 nm and 74 nm, respectively, (a) and (b) show the plane of densities p P , p s , scaled by the 
corresponding bare particle volume v iy i = S y P. The experimental values for average densities are indicated as statepoints (filled circles) with rj s = 0.05 and 
rj P = 0.01, 0.025 and 0.05. The insets illustrate the predictions for the early equilibrium of the experimental mixtures where the nematic (isotropic) regions 
are shown in bright yellow (black), to be compared with the experimental observations after 12 hours (second column of Fig. 1). (c) and (d) show the 
phase diagram in the plane of chemical potentials fip, fi s > scaled by /? = l/(k B T). Sedimentation paths are indicated as red dashed lines. For clarity only the 
paths for rj P = 0.01 (lower) and 0.05 (upper) in (c) and that for rj P = 0.025 in (d) are shown. The square indicates the point of inflection on the binodal. 
The slope of the binodal in (c, d) is related to the slope of the tielines in (a, b) via 21 d/i s /d//pl coex = — Ap P /Ap s ; here Ap z is the density jump of species i. The 
cartoons in (c, d) are to be compared with the experimental observations after 20 months (fourth column of Fig. 1). Also shown is a path for ^ P = 4 mm 
and rj P = 0.05 (green dot-dashed line and upper cartoon). The inset in (d) shows a sedimentation path for a floating nematic in a schematic phase diagram 
with an UN triple point and II demixing that ends in a critical point (circle). 



—oo, strong negative curvature upon increasing fi s and subsequent 
point of inflection. The grand potential of the inhomogeneous sys- 
tem can be decomposed 20 ' 21 as Q = U — TS + a^2 i §Qdzp i (z)(l) exti 
(z) — ifaNh where U is the internal energy, S is the entropy, a 
is the transverse system area, <t>^i{z) = mgz is the external potential, 
z is the height coordinate, H is the sample height, and Nj = 
a§Qdzpj(z) is the total number of particles of species i. Hence 
Q = U — TS — a^2i^Qdzp i (z)iLi i (z), where we combined the bulk 
chemical potential and the gravitational energy into a local chemical 
potential, 

li i (z)=n i -k B Tz/£ i . (1) 

As a b we trust the local density approximation 20 , and obtain the 
density profiles simply via p f (z) = p e ^{fi s {z),fi P {z)). When taking z 
as a parameter, Eq. (1) represents a line with slope £p/£ s in the p P , p s - 
plane. Its endpoints are determined by H (the position of the upper 
meniscus of the dispersion) and the constraints of fixed average 
densities, J ^dzp t (z)/H = rjj Vj. 



As for the small spheres £ s £ Pi the corresponding sedimentation 
paths are almost horizontal, cf. Fig. 2c. Each path crosses only once 
from the nematic to the isotropic side of the binodal. The agreement 
with the experimental findings after 20 months is very good; we took 
partial evaporation of the solvent into account by correspondingly 
decreasing the value of H upon iV z = const. The system with rj P = 
0.05 shows clear compression of the nematic layer, and that with rj P 
= 0.01 possesses a bottom nematic layer, both as experimentally 
observed. 

In the mixture with larger spheres, their larger buoyant mass 
creates significantly steeper sedimentation paths, such that two 
crossings of the binodal can occur, cf. Fig. 2d. The path for rj P = 
0.025 starts in the isotropic region, traverses the nematic, and re- 
enters the isotropic region. This implies the experimentally observed 
floating nematic state. For rj P = 0.01 the nematic region is just 
missed, and the path lies entirely in the isotropic region (not shown). 
However, increasing the value to 0.020, the experimental thin float- 
ing nematic can be very well reproduced (leftmost cartoon in Fig.2d). 
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Figure 3 | Sedimentation profiles obtained from density functional 
theory for mixtures of infinitely thin hard platelets and hard spheres with 
aspect ratio 4.9 (a), and 2.33 (b), representing the effective aspect ratio in 
the experimental mixtures with spheres of 30 nm and 74 nm bare 
diameter, respectively. Shown are density profile for spheres Ps(z)v s 
(dashed black lines) and for platelets p P (z)v P (red dashed lines) as a 
function of height z. The insets show the nematic order parameter S(z) as a 
function of z. Note that the density maximum 5 of platelets in (b) for rj P = 
0.01 is at the location of the experimentally observed thin floating nematic. 
The lower inset in (b) shows density profiles for rj P = 0.020, corresponding 
to the thin floating nematic in the leftmost cartoon in Fig. 2d. 

The agreement for rj P = 0.025 is excellent. For r\ P = 0.05 the point of 
inflection on the binodal facilitates an additional crossing with the 
sedimentation path at high values of fa. The path already starts in the 
nematic region and hence an additional bottom nematic layers 
occurs. As equilibration is expected to be very slow at the dense 
bottom of the corresponding experimental sample, it is not incon- 
ceivable that such a four-layered state would develop after still longer 
waiting times. The slight instability of the lower IN interface in the 
sample with r\ s = 0.025 could be a further hint at such behaviour. 
However, if the thick floating nematic is the equilibrium state for r\ P 
= 0.05, it can easily arise from a slightly steeper sedimentation path, 
e.g. £ P = 4 mm, see Fig. 2d. The floating nematic is a robust feature, 
independent of whether isotropic-isotropic demixing exists or not 
(see inset to Fig. 2d). We present density profiles, p s (z) and p P (z), and 
nematic order parameter profiles, S(z), in Fig. 3. The floating nematic 
is sandwiched between essentially a bottom dense isotropic sphere 
liquid and a top isotropic platelet liquid. 

Discussion 

That particles can accumulate at well-defined heights under gravity 
(or centrifugation) is well-known 5 . In order to relate our theory, 
which is centered around the concept of a local chemical potential, 
to the "generalized Archimedes principle" by Piazza et al 5 , where the 
chemical potential does not explicitly appear, one can use the ther- 
modynamic identity dYl/dpi = J^jPj^f 1 ] j ^Pv when starting from 
Eq. (SI) in Ref. 5 (Supplementary Material); here IT is the pressure 



and the sum is over all components. Applying their theory to the low- 
density limit of the solute (platelet) in the sea of spheres requires a 
test particle calculation in order to obtain the (orientation-depend- 
ent) pair distribution function of spheres around a suspended 
platelet. Our present calculation avoids such computational burdens 
as it works directly on the level of the free energy. Note that the 
platelet density profile shown in the leftmost panel of Fig. 3b (i.e. 
for the lowest value of platelet density), is close to the true low density 
limit of the solute (platelet), which is the object of study in Ref. 5. This 
density profile is very spread out over the entire system, in striking 
contrast to a floating nematic layer. 

A floating phase has been theoretically predicted in colloid-poly- 
mer mixtures 21 . However, the present system facilitates the experi- 
mental observation as £p/£s can he controlled easily by the size ratio 
of the two species. Our considerations generalize to multi-compon- 
ent mixtures, as sedimentation paths will remain straight lines, cf. Eq. 
(1), which cut through binodal (hyper) surfaces. They also apply to 
atomic systems, where f can be of the order of km, and to geo- 
physical scales, where the strength of gravity depends on distance, 
as is relevant e.g. for the radial composition of the Earth 22 . 
Furthermore our results show that care is needed for the correct 
interpretation of phase coexistence under gravity; in particular the 
number of sedimentation layers can exceed the number of coexisting 
bulk phases allowed by the Gibbs phase rule 23 . 
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